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Malaria parasite strains have emerged to tolerate the therapeutic effects of the prophylactics and drugs 
presently available. This resistance now poses a serious challenge to researchers in the bid to overcome 
malaria parasitic infection. Recent studies have shown that FK520 and its analogs inhibit malaria parasites 
growth by binding to FK506 binding proteins (FKBPs) of the parasites. Structure based drug screening 
efforts based on three-dimensional structural information of FKBPs from Plasmodium falciparum led us to 
identify new chemical entities that bind to the parasite FKBP35 and inhibit its growth. Our experimental 
results verify that this novel compound (D44) modulate the PPIase activity of Plasmodium FKBP35 and 
demonstrate the stage-specific growth inhibition of Plasmodium falciparum strains. Here, we present the 
X-ray crystallographic structures of FK506 binding domains (FKBDs) of P/FKBPSS and PvFKBP35 in 
complex with the newly identified inhibitor providing molecular insights into its mode of action. 

With an estimated 800,000 deaths annually, malaria is one of the most prevalent and dangerous infec- 
tious diseases currently afflicting the tropical and sub-tropical regions of India, Africa, Southeast Asia 
and South America''^. Among the five Plasmodium species, P. falciparum and P. vivax, are the most 
virulent and chronic forms that cause malaria infection in humans triggering the lethal cerebral malaria and 
death. In recent years, these parasites have developed resistance to current therapeutics, of particular concern are 
recent reports showing the first appearance of artemisinin resistance in P. falciparum in South East Asia"". This 
poses a serious challenge to the malaria research community, sparking the need for novel protein targets and 
molecules against these drug- resistant forms of malaria. FK506, a FDA-approved drug for organ transplant, 
shows antimalaria efficacy'''^, but its use is contraindicated owing to its potent immunosuppressive actions. 
Subsequently, in search for FK506 binding proteins (FKBPs) in the Plasmodium species, the piggyback approach 
led to the identification of the Plasmodium FKBP35 in both P. falciparum and P. vivax (hereafter referred to as 
P/FKBP35 and PvFKBP35, respectively)*''". 

P/FKBP35 and PvFKBP35 belong to the immunophilin family, more commonly known as FK506 binding 
proteins (FKBPs). Aside from their immunosuppressive properties, which are mediated via interaction with 
calcineurin, members of the FKBP family show a canonical peptidyl-prolyl cis-trans isomerase (PPIase) or 
rotamase activity and regulate various physiological functions, including protein folding, stability'-" histone 
chaperone activity'^'", receptor signaling""", protein trafficking'^"*, transcription''^-''", calcium homeostasis"''^', 
spermatogenesis^"', and neuroprotective and neurotrophic activities^''. These physiological functions emphasize 
the role of FKBPs as helper proteins in assisting and mediating cell responses to physiological signals and effectors 
and therefore emerge as potential molecular targets for pharmacological intervention. However, proof of concept 
study showing the importance of Plasmodium FKBP35 for parasite growth is still elusive. Chemical biology 
approaches have been exploited to verify the role of this protein by employing its high affinity ligand (FK506). Our 
previous microarray study using FK506 and Cyclosporine as chemical probes has revealed strikingly similar 
global transcription profiles of P. falciparum^*. This study points out the role of FK506 binding proteins and 
calcineurin pathway in mediating the parasite growth. Studies from Monaghan et al, have highlighted that the 
inhibition of either PPIase activity or chaperone activity of FK506 binding proteins without calcineurin inhibition 
are sufficient to inhibit the gro^wth of Plasmodium falciparum 3D7 culture^. More importantly, studies from 
Braun et al, have shown that differences in cellular levels of parasite FKBP35 (100 nM) and the host FKBP12 
proteins (5 |iM) offer a unique selectivity window to target parasite FKBP proteins in preference to host proteins'". 
Taken together these studies suggest that Plasmodium FKBP35 could be a viable drug target for controlling 
malaria. 

Of Plasmodium FKBPs, molecular characteristics of P/FKBP35 and PvFKBP35 have been characterized and 
their catalytic and ligand-binding properties are homologous to other FKBP family members, such as FKBPS 1 
and FKBP52^'^''. FK506 binding domain of Plasmodium FKBPs share sequence and structural similarities with 
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other canonical FKBP family members (Supplementary Fig. SI). 
Previously, we have determined the three-dimensional crystal struc- 
tures of P/PKBD35 and PvFKBD35 in complex with FK506, which 
offered insights into the molecular interaction between the FK506 
and Plasmodium FKBD35 proteins^'''^'''^^. In our attempts to discover 
new small molecule inhibitors that could specifically inhibit the para- 
site PPIase activity, we instigated a structure-based in-silico screen- 
ing of commercially available library. Based on our in silico, 
crystallographic and NMR studies, coupled with PPIase and malaria 
growth inhibition assays, here we demonstrate that a purine-like 
derivative devoid of immunosuppressive activity targets Plasmo- 
dium FKBP35, inhibits the parasite growth, and possesses potential 
to be developed as novel therapeutics to combat the malaria. 

Results 

Identification of N-(2-Ethyl-phenyl)-2-(3H-imidazao [4, 5-b] 
pyridin-2-yl-siilfanyl)-acetamide as a novel inhibitor of Plasmo- 
dium FKBPs. Plasmodium FKBD35-FK506 crystal complex struc- 
tures reveal a structural fold, ligand-binding mode similar to those of 
canonical FKBPs^"'^'^. Despite these similarities, the (35-136 loop (H87, 
191) residues in human FKBP12 are varied to Cysl06/105 and 
Serl09/108 residues in Plasmodium species. These residues that 
form the site S2 (Fig. Ic), accommodates the pyranose methyl 
group of FK506, which is further away from the site SI (Fig. lb) 
while the Trp78 residue (site SI) forms the base for the binding of 
pipecolyl moiety of FK506. These varied Cysl06/105 and Serl09/108 



residues are exploited for achieving selectivity towards Plasmodium 
FKBP35. In our structure based pharmacophore modeling, LUDI 
interaction sites (Fig. la) together with excluded volumes (serve to 
provide the shape or volume of protein active site and to prevent 
clashes with protein atoms) highlight the explicit hotspots at the 
active site that are exploited for structure based pharmacophore 
development. These interaction sites are clustered hierarchi- 
cally based on the feature type to minimize the number of features 
that are amenable for pharmacophore development. Two features 
namely - a hydrophobic anchoring Trp78 and an acceptor at Ile75 of 
site SI and similarly three features - an acceptor at TyrlOl, a donor at 
Asp56 and a hydrophobic feature at vicinity of Phe55, Cysl06, and 
Serl09 of site S2, together with exclusion volumes (152) are included 
in model development. Upon screening our in-house ADME 
property filtered 3D virtual database of ChemDiv commercial 
library with this five feature structure based pharmacophore model 
(Fig. Id) has enabled a focused library. Together with molecular 
docking, scoring and consensus results (Supplementary Fig. S2), 
the library were further prioritized based on their interactions at 
above mentioned key residues. One of the top ranking ligand N- 
(2-Ethyl-phenyl)-2-(3H-imidazao [4, 5-b] pyridin-2-ylsulfanyl)- 
acetamide (referred to as D44), which contains a purine-like ring 
system linked to a phenyl ring via thio-acetamide linkage (Fig. 2a), 
selectively inhibits the Plasmodium FKBD35 PPIase activity with an 
IC50 of 132 ± 1.09 nM (Fig. 2b) and 125 ± 2.9 nM for falciparum 
and vivax species respectively. D44 did not affect the enzymatic 




Figure 1 | Structure-based pliarmacophore model, (a) Ludi program identifies discrete acceptor (red rods (OH)), donor (purple rods (NH) and 
hydrophobic (black spheres) vector sites with which a ligand can form efficient interaction with the protein; (b& c) Hierarchical clustering of acceptor, 
donor, and hydrophobic vector site points enables one to select features targeting select residues as depicted at site SI (b) and at site S2 (c). (d) Five 
features were selected targeting specific residues such as F55, D56, 175, W78, and YlOl residues that were included in the structure or receptor based 
pharmacophore (Exclusion volumes are not shown in figure for clarity). 
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activity of either human FKBP12 or FKBP52 even at W \iM 
concentration (Supplementary Fig. S3). The canonical FKBP ligand 
FK506 not only inhibits PPIase activity but also inhibits the 
phosphatase activity of calcineurin''", by forming a ternary complex 
with calcineurin, leading to suppression of immune responses. To 
ascertain that D44 does not interact with calcineurin upon binding to 
FKBP, calcineurin phosphatase assay is performed and our results 
demonstrate that D44 weakly binds to human FKBP 12 but lacks the 
ability to interact with calcineurin and inhibit its phosphatase activity 
(Fig. 2c). 

D44 interferes with P. falciparum intraerythrocytic growth. To 

monitor the effects of D44 on P. falciparum growth, highly syn- 
chronized ring stage (~5 hours post invasion) 3D7 parasites were 
allowed to grow in the presence of D44 inhibitor at concentration 
ranging from 1 nM to 1 |xM for an entire intraerythrocytic deve- 
lopmental cycle (IDC) from the ring stage through egress and 
reinvasion, thus providing readout for assessing the inhibitory 
effects on any part of the erythrocytic cycle. D44 inhibited the 
growth of P. falciparum 3D7 in a dose dependent manner with an 
IC50 of 234.5 ± 1 .32 nM. To map the stage specific inhibitory activity 
of D44, P. falciparum 3D7 parasites growth were closely monitored at 
8 h time interval for an entire IDC in the presence of D44 and FK506 
at ~1 |J,M and 250 nM respectively (corresponding to ~4X IC50). 
Parasites morphology at each time point was monitored on Giemsa- 
stained smears and parasitemia was quantified via flow cytometry. 
From the ring to early schizont maturation, there is no change in 
either parasite morphology (Fig. 3a) or parasitemia (Fig. 3b) between 
the D44- and FK506-treated and control parasites. However dur- 
ing subsequent development, the untreated control parasites 



successively mature to late schizont and ring stages (Fig. 3a) 
resulting to an increased parasitemia (Fig. 3b) whereas both D44- 
and FK506-treated parasites remained arrested at the schizont stage 
(Fig. 3a) and subsequently degraded leading to a decreased para- 
sitemia, suggesting that D44 interferes with the schizont stage 
development similar to FK506^''. 

To validate this finding, we closely monitored the maturation of 
early P. falciparum 3D7 schizont stage parasites (—36-40 hours post 
invasion) in presence of D44 or FK506 and in absence of treatment 
(control). D44 and FK506 treated parasites show significant mor- 
phological differences to untreated parasites (Fig. 3c) only after the 
first 4 hours of schizont maturation. The untreated parasites pro- 
gresses into their subsequent stages (formation of late schizont and 
ring stages) while the D44 and FK506 treated parasites remain 
arrested at the late schizont stage (Fig. 3c) with further degradation 
of D44-treated parasites (Supplementary Fig. S4). This further cor- 
roborates the finding that D44 directly interferes with schizont para- 
sites maturation, probably through a mechanism different from 
FK506 inhibition. 

Furthermore, investigation of the effects of D44 (1 \xM) on an 
additional P. falciparum laboratory- adapted strain (Palo Alto) and 
a field isolate (P/D703) collected from a Senegalese malaria-infected 
patient and adapted to in vitro culture shows that D44 inhibited the 
growth of both strain in a pattern similar to the 3D7 strain (Fig. 3d). 
Since the compounds were dissolved in DMSO, we determined the 
effect of DMSO on parasite growth to discriminate DMSO to com- 
pounds specific inhibition and showed no significant DMSO inhibi- 
tion at 5% dilution (Fig. 3d), suggesting that the observed parasites 
growth inhibition was D44 specific. 
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Figure 2 | (a) Chemical structure of ligand D44. (b) Sigmoid dose response curve showing the inhibition of P/FKBD35 enzymatic activity with the increase 
in concentration of D44. Data is represented as duplicate of two independent experiments ± SEM. (c) Effect of D44 on phosphatase activity of 
calcineurin. The inhibitory effect of FK506, D44 (5 ^M) FKBP12 complexes on the phosphatase activity of calcineurin (CaN) was examined. Calcineurin 
alone serve as a control. Responses were quantified in terms of phosphate released, (d) Dose dependent inhibition of parasitemia of Plasmodium 
falciparum 3D7 cell line. Error bars represent the SEM of two independent experiments in duplicate. 
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Figure 3 | Effects of D44 treatment on the growth and morphology of different Plasmodium falciparum lines (a-b) Morphological characteristics (a) and 
parasitemia quantification (b) of P. falciparum 3D7 parasites growth at 8 h time interval following treatment with either D44 or FK506 showing 
interference of both inhibitors with the schizont stage development compared to control untreated parasites. Data shown in graph (b) are average of 2 
independent experiments performed in triplicate each, (c) Timing and morphological characteristics of the schizont stages inhibition upon incubation 
with D44 and FK506. During the first 4 h, no significant morphological differences were observed between the treated and control parasites. During the 
subsequent development, the control parasites progress to the next generation (formation of mature schizont and ring) while the treated cells remain 
arrested at the late schizont stage, (d) Parasitemia determination of laboratory-strain (P/Palo Alto) and field isolate (P/0703) after a complete IDC 
(48 hrs) in presence of D44 or controls DMSO-treated and untreated parasites. D44 significantly inhibits growth of both parasites strains similar to 3D7 
line. 



Crystallographic structures of Plasmodium FKBD35 in complex 
with D44. To gain better structural insights into the mode of action 
of D44, X-ray crystallographic structures of P/PKBD35 and 
PvFKBD35 in complex with D44 were determined to a resolution 
of 2.75 A and 1.73 A, respectively. Refined crystallographic structure 
revealed a topology consisting of six P-stranded sheets and a short i- 
helix similar to the previous Plasmodium FKBD35-FK506 structures 
from P. falciparum (PDB ID aVNl)"* and P. vivax (PDB ID 3IHZ)'''. 
The P/PKBD35-D44 complex consists of three monomers in the 
asymmetric unit (Fig. 4a) with chains A and B forming a covalent 
dimer via an intermolecular disulphide bond between Cysl06 from 
both the chains. Similarly the third monomer (chain C) also forms an 
identical dimer with the neighboring chain, which is generated by 
crystal packing. While PvFKBD35-D44 complex consisted of one 
monomer in the asymmetric unit (Fig. 4b) forming a dimer with 
the neighboring chain, generated by crystal packing, with the help 
of a disulphide bond through CyslOS. The Pf and PvFKBD35 
structures superpose well with a backbone root-mean-square devia- 
tion (RMSD) of 0.314 A. 

Discussion 

Our ligand design was mainly aimed to target the two highly con- 
served Cysl06/105 and Serl09/108 residues within the P. falciparum 
and P. vivax species to achieve selectivity towards Plasmodium 
FKBP35 (Supplementary Fig. 81). Virtual screening with structure 
or receptor based pharmacophore model with exclusion volumes are 




Figure 4 | The crystal structures of Plasmodium FKBD35 in complex 
with the ligand D44. The cartoon representation of the asymmetric units 
of P/FKBD35 (a) and PvFKBD35 (b) in complex with ligand D44 (in stick 
mode), showing the three chains and a single chain respectively. 2Fo-Fc 
difference Fourier map contoured at 1 a cut-off of ligand D44 in complex 
with P^KBD35 (c) and PvFKBD35 (d), in stereo mode. Figures were 
generated using the program PyMOL''^. 
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shown to map hits that are complementary to the active site and 
improve the hit rate by filtering the false positives". Ligands which 
are predicted to have interaction at both sites SI and S2 and in 
consensus with most of the scoring functions are selected for experi- 
mental screening. Inhibition of substrate activity of the PPIase pro- 
vides a convenient assay to evaluate whether the ligands bind and 
inhibit the enzymatic activity of the protein. D44 potently inhibits 
the PPIase activity of P.falciparum/vivax FKBP35 at sub-micromolar 
concentration (132 nM/125 nM) selectively and show —100 fold 
selectivity towards parasite FKBPs over human FKBPs. D44 not only 
inhibits the in vitro growth of Plasmodium falciparum 3D7 strain 
parasite at 234.5 ± 1.32 nM (Fig. 2d) but also the growth of Palo Alto 
strain and field isolate of P/D703 (Essai 15/05/2013, collected from a 
Senegalese malaria infected patient) by 50% at 1 [iM concentration. 
These data suggests that D44 can inhibit the PPIase activity of 
FKBP35 at sub-micromolar range and hinder the growth of different 
strains of Plasmodium falciparum at micromolar range. Further lead 
optimization studies are needed to improve the efficacy of its anti- 
malaria properties. D44 exhibited low cytotoxicity on HeLa and 
HTC116 cancer cell lines and their cell proliferation was inhibited 
albeit at much higher concentrations with an IC50 of 465 ±3.16 [iM 
and 343.6 ± 1.26 |iM respectively, suggesting 1000 fold selectivity 
towards malaria parasite (Supplementary Fig. S5). It is important to 
note that direct evidence linking the inhibition of PPIase activity to 
growth of parasite is lacking, nonetheless many studies on human 
orthologues have shown the importance of FKBPs in regulation of 
cell cycle"", transcription"'^" and histone methylation'^'". Further 
studies using chemical biology approaches using selective Plasmo- 
dium FKBP ligands like D44 are awaited to delineate the role of 
FKBPs or its substrate protein interactions that might be critical to 
the growth and survival of the malaria parasite. 



Comparison to previously determined apo form NMR struc- 
tures'"'^'" indicate a well conserved structural fold with an RMSD of 
2.28 A and 1.37 A for P/FKBD35 and PvFKBD35, respectively 
(Supplementary Fig. S6a-b). As anticipated, the ligand flanking 
loops namely (33 -(54, P4-ck1, and (55 -(56 undergo small structural 
changes leading to an increase in accessible surface area 
(Supplementary Fig. S6a-b). The psi angle of Phe55 undergoes a 
change from 124° (apo form) to —21° (D44) which is 6°/2.6° more 
than the FK506 complex in Plasmodium falciparum/ vivax structures 
respectively. These changes could well be attributed to the non- 
bonded contacts made by D44 with Phe55/54, unlike the hydro- 
gen-bonded contacts of FK506, which gives the residues more 
flexibility. The electron density for D44 (Fig. 4c, d) was traced 
unequivocally in the electron density difference map in the hydro- 
phobic binding site. One ligand per monomer was observed, sitting 
in the dimer interface facing each other separated by a distance of 
3.6-4.2 A. Crystallographic structures ofP/FKBD35 andPvFKBD35 
in complex with D44 showed that the purine like ring of D44 makes 
critical hydrogen bonding contacts with the catalytic residues Asp56/ 
55 and TyrlOl/100 (Fig. 5a, b) and the thio-acetamide linker of D44 
through its carbonyl group interacts with Ile75/74, locking the ligand 
in the active site pocket (Fig. 5a, b & c) (Table SI). Further non- 
bonded contacts (Fig. 5a, b) with residues like Phe55/54, Phe65/64, 
Trp78/77, and Ilel 10/109 stabilize the ligand-bound complexes 
(Table SI). 

Our 2D 'H,'''N heteronuclear single quantum correlation (HSQC) 
nuclear magnetic spectroscopy (NMR) results confirmed that most 
of the active site residues including the Phe65, Trp78 residues show 
chemical shift perturbation upon binding with D44 (Supplementary 
Fig. S7). Superimposition FKBD35-D44 structures to apo form struc- 
tures*"""" not only reveal a translation of —2-6.8 A (Fig. 5d) on 
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Figure 5 | Comparison of binding site interactions, (a & b) The stick representation of residues making hydrogen bonds (black) and the non-bonded 
contacts (brown) with ligand D44 in P/PKBD35 (a) and PvFKBD35 (b). Ser 109/108, one of the highly conserved residues is also shown for clarity; (c) 
Electrostatic surface representation of P/PKBD35 in complex with D44 (in stick mode), showing the docking of the ligand into the hydrophobic active site 
cleft; (d & e) Superposition of interacting residues of apo (P/ - light orange; Pv - orange) FKBD35 with their corresponding FK506 (P/- light green; Pv - 
green) and D44 (Pf- light blue; Pv - dark blue) complexes. The D44 of P/PKBD35 alone is also shown for clarity. Large deviations in hydrogen bonded 
residues can be observed for PjPKBDSS compared to its P.vivax counterpart (d). While the residues involved in non-bonded contacts (e) clearly show a 
movement away from the active site as the ligand gets bigger, compared to the apo form; (f) Specificity of ligand D44 towards the highly conserved Cys 
106/105 and Ser 109/108 can be observed when superposed with the corresponding FK506 structures. D44 (P/- light green; Pv- dark green) clearly moves 
closer to Cys 109/108 compared to FK506 (Pf- light blue; Pv- cyan). The distances are shown for the PvFKBD35, with D44 at 3.76 and 3.74 A from Cys 
105 and Ser 108 respectively, compared to 4.04 and 8.05 A in FK506. In addition, the hydroxyl group of Serl08 flips towards the ligand (b & e) in 
PvFKBD35 structure enhancing the ligand specificity. It could also be noticed that the thio-acetamide linker in D44 compensates for the pipecolinyl 
moiety (shown within the black open circle) of FK506, a unique feature pertaining to D44. In figure 5 d, e & f the residues are labeled according to their 
species, namely Pf/Pv. The figures were generated using the program PyMOL"^. 
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residues involved in hydrogen bonded interactions, but also notice- 
able changes in the local conformations of Phe65/64 and Trp78/77 
side chains (Fig. 5e). Quantitative comparison of ligand bound struc- 
tures of Plasmodium FKBDs indicate that the Phe65/64 moves by 
-4/2.8 A and -5.5/3.5 A in complex with D44 and FK506, respect- 
ively, while Trp78/77 sways by ~6/5.7 A and -7.5/5.9 A (Fig. 5e). In 
addition, directionality can be conferred as these two residues appro- 
ximate towards the thio-acetamide due to the strong hydrophobic 
interactions of sulphur with the hydrophobic residues (Fig. 5e). In 
P/FKBD35, D44 is much closer (Fig. 5f) to the nearest atom from the 
purine like ring located at a distance of 3.6 A and 4.3 A from Cysl06 
and Serl09 respectively, compared to 3.9 A and 5.4 A in the corres- 
ponding FK506 complex. Similarly in PvFKBD35, CyslOS and 
SerlOS are at a distance of 3.76 A and 3.74 A from D44, compared 
to 4.04 and 8.05 A in its FK506 (Fig. 5f) complex. The hydroxyl group 
of Serl08 in PvFKBD35 flips towards purine like ring of D44 (Fig. 5b, 
f) due to close van der Waal contacts with this residue which is well 
conserved among the Plasmodium species, thereby enhancing its 
selectivity towards Plasmodium FKBPs. Structure activity relation- 
ship studies on FKBP12 inhibitors^^'^^ have shown the presence of 
pipecolinyl or proline moiety being paramount for the binding and 
inhibition of enzymatic activity. Here, we present a unique chemo- 
type which possesses a thio-acetamide compensating the lack of 
bulky pipecolinyl group (Fig. 5f) and through its hydrophobic inter- 
actions induce local conformational changes (Fig. 5e, f) at the active 
site. Overall, though D44 (MW of 312) is smaller in size compared to 
FK506 (MW of 804), it interacts with specific residues as predicted, 
confirming that these interactions between the protein and ligand are 
critical for the inhibition of the enzyme. 

Taken together, we have identified a novel small molecule Plasmo- 
dium FKBP35 antagonist and demonstrated its biological activity in 
P. falciparum. The elucidation of crystal structures of Plasmodium 
FKBD35 in complex with D44 further substantiates the efficiency of 
ligand to target not only P. falciparum but also P. vivax, a much 
neglected target whose incidence is re-emerging with drug resistant 
forms and gaining awareness^^^'^^. Thus, the information presented in 
the current study could be used to aid designing selective small 
molecule inhibitors that bind to the parasite FKBP and potentially 
antagonize their function in the parasite life cycle. 

Methods 

Structure-based pharmacophore. P/FKBD35-FK506 X-ray crystal structure (PDB 
ID 2VN1) was prepared by correcting bond orders, missing atoms and energy 
minimized using conjugate gradient force field. Ludi interaction map^^"^^ was 
generated around 10 A of the active site residues using three features - acceptor, donor 
and lipophilic features highlighting the probable interaction sites for the ligands with 
the protein atoms. Features anchoring residues like Asp56, Trp78, TyrlOl, and 
Cysl06 are included in structure based pharmacophore model. To avoid clashes of 
ligand with the protein atoms, a total of 1 52 exclusion volumes are incorporated into 
the structure based pharmacophore modeP^. Pharmacophore screening on ChemDiv 
library resulted in a focused library of 13,000 molecules. This library is further 
processed by applying ADME filters which refmed the library to —2600 molecules. 

Virtual screening and consensus scoring. An advanced molecular docking program 
GOLD v4.0^^ was employed in docking the focused library into the active site of 
P/FKBP35. Docking parameters were optimized to simulate the crystal ligand 
conformation and the same parameters were used for screening the focused library. 
Thus obtained docking solutions were not only scored with default gold scoring 
function, but also rescored with Ligscore'^ PLP*", ]AIN^\ LUDP^ and PMF^-*' 
scoring functions in Accelrys suite. A consensus scoring scheme**"*^ was applied to 
remove any bias or inadequacies of using one particular scoring function. Upon visual 
inspection of 500 top scoring molecules, compounds were prioritized for contact to 
key catalytic residues like Asp56, TyrlOl, and Ile75 and accordingly 46 compounds 
were purchased and tested experimentally {Supplementary Fig. S2). 

PPIase assay. PPIase activity was assessed in the presence and absence of inhibitors 
by the method described earlier^" using Succ-Ala-Leu-Pro-Phe-p-nitroaniline as a 
substrate (Peptide Institute, Osaka, Japan). All the reagents were pre-equilibrated at 
O^C The final reaction mixture {1 mL) contained 100 nM of P/FKBD35 or human 
FKBP12/FKBP52 (100 nM), 100 |iM of substrate and 300 |ig of a-chymotrypsin 
(Sigma, St. Louis, MO) and assay buffer (50 mM HEPES, 100 mM NaCl, pH 8.0 at 
0°C). Reaction was performed at 4"C and the changes in the absorbance were 



monitored at 390 nm for 3 min on a Shimadzu UV-1800 UV-vis spectrophotometer. 
IC50 was calculated from dose response curve in duplicate using GraphPad Software^'. 

Calcineurin assay. Calcineurin assay was carried out using recombinant human 
calcineurin as described previously^^"^^, with some modifications. Reactions (50 |j.l) 
contained 50 mM Tris-HCl (pH 7.5), 100 mM NaCl, 6 mM MgCla, 0.5 mM CaCU, 
0.1 mg/ml BSA, 0.5 mM DTT, 0.025% NP-40, 0.25 ^iM calmodulin, HsFKBP12, 5 
|iM FK506 and D44, 0.5 |iM RII phosphopeptide (a specific calcineurin substrate). 
Reaction mixtures were incubated for 10-30 min at 30 C before the addition of the 
substrate peptide to initiate dephosphorylation at 30 'C. The reaction was terminated 
by addition of Malachite Green. The reaction was monitored by reading the 
absorbance at 620 nm. 

Growth inhibition assays. P. falciparum laboratory strains (3D7 and Palo Alto) and 
field isolate {P/D703) were used in this study to investigate the effects of D44 on 
parasites growth. The P/0703field isolate was isolated from a Senegalese malaria- 
infected patient and adapted to in vitro in the laboratory at the Pasteur Institute of 
Dakar. Parasites (5% parasitemia, 2% hematocrit) at the mid ring stage (—6-10 hours 
post-invasion, hpi) were incubated in triplicate in a 96 well plate with 250 nM of 
FK506 or 1 |j,M of D44 or in absence of treatment (control). Parasites maturation was 
monitored every 8 h for an entire Intraerythrocytic Developmental Cycle (IDC) by 
Giemsa staining of fixed smears for comparison of parasite morphology at each time 
point. Parasitemia was quantified by flow cytometry based assay using either 
Ethidium Bromide (1/1000) (for 3D7) or Sybr Green I (1/1000) (for Palo Alto and 
P/D703) DNA stain dyes. Briefly, parasites were incubated with selected either dye in 
PBS for Ihour at room temperature, washed 3 times with PBS then resuspended in 
PBS. Wells content were separately transferred to a FACS tube and parasitemia was 
quantified through a FACS Calibur with an acquisition of 10.000 events/sample. 

Monitoring of stage specific effects of D44 and FK506 was carried out by incubating 
separately highly synchronized ring, trophozoite or schizont stages of P. falciparum 
3D7 parasites with FK506 (250 |iM), D4 (1 |iM) or without inhibitor (control) for 
24 h, duration sufficient to allow transition to the next stage. After 8 h of incubation, 
drugs were removed, parasites were washed twice with PBS, then fresh growth 
medium was added into the culture and parasites were allowed to mature for an 
additional 16 h. Parasite morphology and parasitemia were determined as described 
above. 

Cell proUferation assay. HeLa, HTC116 cell lines were maintained in Dulbecco's 
Modified Eagle Medium (DMEM) containing 10% fetal bovine serum, 1% penicillin/ 
streptomycin and maintained at 37'C in a humidified 5% CO2 chamber. Cell 
proliferation was determined using the CellTiter 96® AQueous One Solution Cell 
Proliferation Assay kit (Promega, Madison, WI, USA). This assay uses a colorimetric 
method to determine the number of viable cells. Briefly, 5X10^ cells/ well were seeded 
in a 96-well plate. The volume of the media in each well was 100 ]A. Cells were 
maintained in a humidified, 5% CO2 atmosphere tOl the monolayers reached 50-60% 
confluency (24-36 h). Cells were then incubated with D44 for another 72 h. At the 
end of this period, 20 )j.l/weU of the reagent was added to the cells, incubated for 
another 1-4 h and absorbance was measured at 490 nm in a 96-weU plate reader. The 
values shown are the mean ± SEM of at least three independent experiments 
performed in tripUcates. IC50 were calculated based on A490 (as % control of 
untreated cells). 

Crystallization, X-ray data collection, structure determination and refinement. 

P/FKBD35 and PvFKBD35 were expressed and purified as described previously^^ and 
was used for crystallization at a concentration of 15 mg/ml. Initially, a series of 
ligands were screened for crystallization as a cocktail containing two or three ligands; 
the cocktaU approach proved screening for initial crystal hits efficient and 
comprehensive^^. The cocktail of ligands were incubated overnight with the protein at 
a molar ratio of 1 : 1. This protein- cocktail complex was screened using an 
ammonium sulphate-buffer grid. After 2 weeks, rod shaped crystals appeared in 
3.5 M ammonium sulphate and 0.1 M citrate buffer, pH 5.5 with P/FKBD35-D44 
containing cocktail. Further optimization with D44 alone yielded bigger crystals in 
3.5 M ammonium sulphate and 0.1 M citrate buffer, pH 5.75. Similarly, the 
PvFKBD35-D44 crystals grew in 3.0 M ammonium sulphate and 0.1 M bicine, pH 
9.0. Both crystals were cryo- protected with 20% glycerol added to their corresponding 
mother liquors and flash-frozen at 1 00 K during data collection. Data for P/FKBD35- 
D44 from a single crystal up to 2.75 A were collected at 100 Konbeamline 13B1 at the 
National Synchrotron Radiation Research Center (Hsinchu, Taiwan) using ADSC- 
Quantum 315 detector at a wavelength of 0.9642A. While data up to 1.73 A was 
collected from a single crystal of PvFKBD35-D44 using Rigaku R-Axis IV image plate 
mounted on a Rigaku Micromax-007 HF X-ray generator with CuK" radiation. The 
diffraction data was indexed, integrated and scaled using the HKL2000^^ suite of 
programs and d*TREK^^. The PfFKBD35-D44 crystal belonged to the tetragonal 
space group P432i2, similar to the previously solved P/FKBD35-FK506 complex^^ 
(PDB ID 2VN1). The Mathews co-efficient indicated the presence of three monomers 
in the asymmetric unit, due to an enlarged unit cell. Therefore, a molecular 
replacement search was employed using the program PHASER^", with a monomer 
unit of the P/FKBD35-FK506 (PDB ID 2VN1). PHASER identified the three 
monomers in the asymmetric unit. A similar protocol was adapted for PvFKBD35- 
D44 data which belonged to the trigonal space group P3221, with one molecule in the 
asymmetric unit located using the program PHASER with a monomer unit of 
PvFKBD35-FK506 (PDB ID 3IHZ)^^ as the search model. In both structures, after 
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initial refinement using REFMAC^\ the electron density was observed for D44 in the 
Fo-Fc difference Fourier map at 3 a cut-off. Further cycles of map fitting using the 
program COOT^^, followed by refinement, resulted in the final structure with good 
geometry and stereochemistry. A few missing residues were observed at the terminal 
ends of the protein. The program MOLPROBITY^^ showed none of the residues 
outside the Ramachandran plot^*. The data collection and refinement statistics are 
given in Table S2. All structural data, including structure factors where applicable, 
have been deposited into the Protein Data Bank with PDB codes 4J4N & 4J40. 

Chemical shift perturbation analysis by NMR spectroscopy. To study the ligand- 
protein interactions, chemical shift perturbations to ^^N- labeled Plasmodium 
FKBD35 were monitored on 2D HSQC spectra recorded at 298K on a Bruker Avance 
700. For this, 0.1 mM protein samples in 20 mM sodium phosphate, pH6.5 were used 
and chemical shift perturbations were analyzed by adding an excess amount of the 
ligand D44. 
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